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[ij  We  use  recently  acquired  geomagnetic  archival  data  to  extend  our  long-term 
reconstruction  of  the  heliospheric  magnetic  field  (HMF)  strength.  The  1835-2009  HMF 
series  is  based  on  an  updated  and  substantiated  IntcrDiumal  Variability  (IDlr)  series  from 
1872  onwards  and  on  Bartels’  extension,  by  proxy,  of  his  //-series  from  1835  to  1871.  The 
new  IDV series,  termed  1DV09,  has  excellent  agreement  (R2  =  0.98;  RMS  =  0.3  nT)  with 
the  earlier  1DV05  series,  and  also  with  the  negative  component  of  Love’s  extended  (to 
1905)  Dst  scries  (R2  =  0.91).  Of  greatest  importance  to  the  community,  in  an  area  of 
research  that  has  been  contentious,  comparison  of  the  extended  HMF  series  w  ith  other 
recent  reconstructions  of  solar  wind  B  for  the  last  -100  years  yields  a  strong  consensus 
between  series  based  on  geomagnetic  data.  Differences  exist  from  -1900-1910  but  they 
are  far  smaller  than  the  previous  disagreement  for  this  key  interval  of  low  solar  wind  B 
values  which  closely  resembles  current  solar  activity.  Equally  encouraging,  a  discrepancy 
with  an  HMF  reconstruction  based  on  l0Bc  data  for  the  first  half  of  the  20th  century 
has  largely  been  removed  by  a  revised  mBc-based  reconstruction  published  after  we 
submitted  this  paper,  although  a  remaining  discrepancy  for  the  years  -1885-1905  will 
need  to  be  resolved. 
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I.  Introduction 

[:]  Svalgaard  and  Cliver  [2005]  introduced  the  Inter- 
Diurnal  Variability  ( IDV)  index  for  a  given  geomagnetic 
observatory  (“station")  as  the  average  difference  without 
regard  to  sign,  from  one  day  to  the  next,  between  hourly 
mean  values  of  the  Horizontal  Component,  //,  and  measured 
one  hour  after  midnight.  The  average  should  be  taken  over  a 
suitably  long  interval  of  time,  such  as  one  year,  to  eliminate 
various  seasonal  complications. 

[3]  IDV  has  the  useful  property’  of  being  independent  of 
solar  wind  speed  and  is  highly  correlated  with  the  near- 
Earth  heliospheric  magnetic  field  (HMF)  strength  B.  Thus 
once  IDV  is  determined,  solar  wind  B  is  known  as  well. 
Svalgaard  and  Cliver  [2005]  used  IDV  augmented  with 
Bartels '  [1932]  //-measure  to  reconstruct  the  HMF  strength 
for  the  years  1872-2004. 

[4]  1  lere  we  report  on  an  extension  of  the  IDV  index  for  a 
longer  time  interval  ( 1 835  2009),  using  many  more  stations. 
The  inclusion  of  more  data  is  particularly  important  for  the 
years  from  1890  to  1909  for  which  the  initial  version  of  the 
index  (1DV05)  was  based  on  observations  from  only  one 
station  before  1901  and  four  more  stations  from  1903.  An 
important  aspect  of  1DV09  is  that  it  includes  recent  years 
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with  index  values  at  the  same  level  as  the  very  low  values  in 
1901-1902,  thus  allowing  the  correlation  between  IDV  and 
the  magnitude  of  the  near  Earth  I  IMF  to  be  extended  to  such 
low  values  without  extrapolation.  With  this  correlation,  we 
infer  HMF  B  for  years  prior  to  the  space  age  and  compare 
our  B  values  with  those  obtained  by  other  investigators 
using  geomagnetic  or  cosmic  ray  data. 

2.  Analysis 

2.1.  Derivation  of  IDV 09 

[5]  Our  determination  of  1DV09  is  essentially  identical  to 
that  of  IDV05  except  for  the  inclusion  of  more  data. 
Svalgaard  and  Cliver  [2005]  emphasized  that  IDV  is  a 
modem  version  of  the  //-'measure  building  on  ideas  of  a 
century  ago  [Moos,  1910].  Kertz  [1958],  May  and  [1980], 
and  Svalgaard  [2005]  suggested  using  only  nighttime  values 
to  avoid  contamination  by  the  regular  diurnal  variation,  SH. 
We  followed  their  lead  but  further  limited  the  time  interval 
to  only  one  hour  following  local  midnight  and  constructed 
the  IDV  for  a  given  station  as  the  unsigned  difference 
between  two  consecutive  days  of  the  average  value  over  the 
interval  of  a  II  component  measured  in  nT.  The  individual 
unsigned  differences  were  then  averaged  over  longer  intervals, 
c.g.,  one  ftill  year  (minimizing  various  geometric  and  sea¬ 
sonal  effects,  e.g„  the  semiannual  non  solar  variation  due  to 
the  tilt  of  the  Earth's  dipole  a  plot  of  27-day  Bartels 
Rotation  values  of  IDV  is  given  by  Svalgaard  [2009]). 

[ft|  Since  2005,  we  have  been  collecting,  digitizing, 
quality  controlling,  and  correcting  (where  needed)  hourly 
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Tabic  1.  Stations  Used  for  IDV09,  Including  Replacement  Sta¬ 
tions  due  to  Relocation  of  Original  Stations8 

Corrected 


Stations 

Geodetic 

Latitude 

Geodetic 

Longitude 

Geomagnetic 
Latitude  2000 

Divisor 

BOX 

58.0 

39.0 

53  9 

1  10 

1  SKb 

55.3 

356.8 

52.9 

1.00 

FKT,  SVD,  ARS 

56.4 

58.6 

52.1 

1  10 

RSV,  BFE 

55.6 

11.7 

52.1 

1  10 

MOS 

55.5 

37.3 

51  3 

1  10 

NVS 

55.0 

829 

50.5 

0.97 

WLH.WNG 

53.7 

9.1 

50.1 

097 

MNK 

54,1 

26.5 

49  9 

0  98 

CLH, FRD 

38.2 

282.6 

49  7 

0  97 

BOU 

40.1 

254  8 

49  2 

0.99 

BAL 

38.8 

264.8 

49  0 

0.99 

DBN,  W  IT 

52.1 

5.2 

48,4 

0.98 

lOu 

524 

13.1 

48  3 

1.00 

POT,  SED,  NGK 

52.1 

12.7 

48.0 

1.00 

\BN,  HAD 

51.0 

355.5 

47.8 

0.99 

BEL 

51.8 

20,8 

47  5 

1  01 

ir  r 

522 

104.5 

47  0 

1  02 

TKT 

41  3 

69  6 

46  5 

1.08 

PET 

53.1 

158  6 

46  3 

1.02 

rx)u 

50.1 

46 

46  0 

1.02 

LVV 

49  9 

23,8 

45.3 

1.04 

PSM,  VU,  CLF 

48.0 

2.3 

43  6 

1.04 

FUR 

48  2 

11  3 

43  4 

1.05 

HRB 

47  9 

18.2 

43  0 

1.06 

THY 

46  9 

17.9 

41  8 

1  08 

YSS 

47.0 

142.7 

39  9 

1  10 

rue 

32.3 

249  2 

39.9 

1  10 

AAA 

43.3 

76  9 

38  4 

1  12 

TFS 

42.1 

447 

37.2 

1  14 

MMB 

43  9 

1442 

36.7 

1  13 

AQU 

42  4 

13  3 

363 

1  13 

BJI.  BMT 

40.3 

1  16.2 

34.2 

1  16 

SFS,  EBR 

40  8 

0.5 

34  2 

1.14 

COI 

402 

351.6 

34.1 

1  15 

LNP,  L7H 

36.1 

103.9 

30  1 

1  20 

VQS.  SJG 

184 

293  9 

29.2 

1  20 

KAK 

36.2 

140  2 

28  9 

1  20 

KNZ 

35  3 

1400 

27.9 

1  21 

HTY 

33.1 

1 39.8 

25.7 

1.23 

SSH 

31.1 

121.2 

244 

1.24 

KNY 

31.4 

130.9 

24.3 

1.24 

HON 

21.3 

202.0 

21.7 

1.26 

GUI 

28.3 

343.6 

15.7 

1.29 

PHU 

21  0 

106.0 

13.7 

1  30 

ABG 

18.6 

72.9 

11  8 

l  31 

KOU 

5.1 

307.4 

108 

1  30 

MBO 

14.4 

343  0 

3  2 

1.31 

ANN 

11.4 

79  7 

3.1 

1  32 

TAM 

22.8 

5.5 

3.1 

1.32 

HUA 

-12.1 

284,7 

2.1 

1  32 

GLIA 

13.6 

144.9 

1.0 

1  32 

IRD 

8.5 

77.0 

0.4 

1  32 

AAE 

9,0 

38.8 

-1.3 

1.32 

BNG 

44 

18.6 

22 

1  32 

asc 

-7.5 

345  6 

-7.9 

1  32 

btv 

-6.2 

1068 

15.8 

1  29 

PPT 

176 

210.4 

-164 

1  29 

vss 

22  4 

3164 

16,5 

1  30 

PIL 

31.7 

296.1 

-18.6 

1.28 

API 

-138.8 

188.2 

-24.6 

1.25 

TAN 

18.9 

47.6 

-29.1 

1.20 

rsu 

-19.2 

17  7 

-30  0 

1  20 

II  UK 

22  9 

27.7 

33  6 

1.17 

CTO.  HER 

34.4 

19.2 

-42.3 

1  09 

WAT,  GNA 

31.8 

1 16.0 

44.4 

1  05 

KX).  (  NB 

35  3 

1494 

-45  8 

1.04 

IRW 

43.3 

19  0 

47  8 

1  02 

Table  I.  (continued) 


Stations 

Geodetic 

Latitude 

Geodetic 

Longitude 

Corrected 
Geomagnetic 
Latitude  2000 

Divisor 

A  MS1* 

-37.8 

77.6 

49.1 

1  00 

AlA 

-65.2 

295.7 

-49.8 

1  20 

AML.  EYR 

-43  4 

172.4 

-50.3 

0  97 

CZT 

-46  4 

51.9 

53  1 

1.10 

aThe  Corrected  Geomagnetic  Latitude  for  the  year  2000  is  given  for 
illustration,  but  the  centroid  of  the  latitudes  for  the  time  of  operation  was 
used  to  estimate  the  Normalization  Constants.  Constants  in  italics  were 
determined  by  an  empirical  fit  to  time  overlapping  stations  Stations  are 
given  by  the  IAGA  abbreviation.  A  list  of  1AGA  designations, 
observatory  names,  and  other  station  details  can  be  found  at  http://www. 
geomag  bgs  ac  ulc  gifs  annual  means  shtml. 

bWeakly  nonlinear  relationships  have  been  used  to  normalize  directly  to 
NGK. 


historical  geomagnetic  data  from  individual  observatories  as 
well  as  from  World  Data  Centers  (there  is,  as  yet,  no 
mechanism  available  to  individual  researchers  for  injeettng 
new  or  corrected  data  into  the  World  Data  Centers  or  into 
various  National  Depositories,  so  we  offer  the  data  to 
interested  researchers  upon  request).  Here  we  use  these 
newly  acquired  data  to  substantiate  the  /DC-index,  which  is 
especially  important  for  the  first  ~30  years  of  the  time  series 
(1872-1902),  during  which  IDV05  was  based  solely  on 
Schmidt's  [1926]  and  Bartels'  [1932]  //  measure  from  1872 
to  1889,  on  Potsdam  observations  from  1890  to  1902,  plus 
Cheltenham  for  1901-1902,  and  Honolulu  for  1902.  In 
contrast,  IDV09  is  based  on  tour  times  as  many  “station 
years”  (143  versus  34)  for  this  31-yr  interval  as  detailed  in 
Data  Set  SI  of  the  auxiliary  material.1  We  update  the  time 
series  by  adding  the  index  values  for  2004-2009.  These 
latter  years  are  significant  because  the  yearly  averages  of  B 
observed  in  2007-2009  are  the  lowest  observed  during  the 
space  age.  They  lie  at  the  lower  endpoint  of  the  correlation 
between  yearly  averages  of  observed  B  and  IDV. 

[7]  Table  1  contains  a  list  of  the  71  stations  (including 
replacement  stations)  used  to  compute  1DV09  (versus  14  for 
IDV05).  A  comprehensive  list  of  the  data  coverage  and  the 
data  values  for  the  individual  stations  used  in  this  study  is 
given  in  Data  Set  SI  in  the  auxiliary  material.  All  raw  data  is 
available  from  the  authors  (LS)  upon  request. 

2.1.1.  Latitude  Normalization 

[h]  For  1DV05,  we  normalized  IDV  values  for  a  given 
station  with  Corrected  Geomagnetic  Latitude,  AY,  to  those  of 
Niemegk  (NGK;  as  Bartels  did  for  the  //-measure)  using 

IDVmvm  inVm/(  1.324  cos"  7(W))  (I) 

[9]  Here  we  have  retained  this  relationship  for  stations 
with  |jV/|  <51°  because  it  still  fits  the  data  for  the  additional 
stations.  At  significantly  higher  latitudes,  the  index  becomes 
strongly  contaminated  by  auroral  zone  activity  [see 
Svalgaard  and  C liver ,  2005,  Figure  2],  and  we  recommended 
not  using  such  stations,  e.g.,  the  long-running  station 
Sodankyla,  SOD  (used  by  Lockwood  et  al.  [2009]).  For 

Auxiliary  materials  are  available  at  ftp  ftp  agu  org  apend  ja 
2009ja0l  5069 
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Normalization  Divisor 


Figure  1.  Adopted  divisors  (blue  circles)  to  normalize  IDV 
to  the  NGK-scale  as  a  function  of  average  corrected  geo¬ 
magnetic  latitude  for  each  station  over  the  time  of  operation. 
For  each  station,  different  color  coded  symbols  show  what 
the  divisor  would  have  been  for  that  station  for  years 
1800  (pink  pluses),  1900  (orange  triangles),  and  2000  (red 
diamonds), 


IDV09,  we  relax  this  restriction  slightly  (by  a  few  degrees 
for  a  few  stations,  indicated  in  Table  I  using  a  constant, 
empirical  normalization  divisor  of  1.1  for  these  stations, 
instead  of  the  divisor  in  equation  (I)).  A  value  larger  than 
-0.95  for  |A/|  >  51°  indicates  some  contamination  by 
auroral  zone  activity.  We  have  not  attempted  to  further 
quantify  the  latitude  dependence  of  the  contamination,  but 
simply  use  an  average  value  for  the  few  stations  slightly 
above  51°.  We  do  this  to  accommodate  changes  in  M 
with  time  which  for  some  stations  can  exceed  several 
degrees  and  to  include  a  few  long-running  stations  just 
above  51°.  (We  expect  only  a  very  weak  influence  in  the 
basic  response  of  the  Ring  Current  (see  section  2.1.5)  to 
the  change  of  the  Earth’s  magnetic  dipole  moment  (as  per 
Glassmeier  et  al  [2004])  over  the  interval  in  question, 
and  so  have  not  attempted  to  correct  for  this.)  Figure  I 
shows  the  adopted  normalization  divisor  as  a  function 
of  W  for  the  71  stations  used  in  the  present  study.  Dif¬ 
ferent  symbols  denote  the  divisor  values  for  the  years 
1800,  1900,  and  2000,  showing  the  sensitivity  of  IDV  to 
changes  in  latitude.  The  normalization  divisor  was  cal¬ 
culated  for  the  centroid  of  the  latitudes  for  the  actual  data 
coverage  for  each  station.  If  we  did  not  normalize,  the 
presence  of  data  gaps  (of  which  there  are  many)  would 
produce  discontinuities  in  the  composite  series. 

2.1.2.  Effect  of  Hourly  Means  Versus  Hourly  V  alues  on 
IDV 

[to]  Farly  magnetometer  data  were  taken  (and/or  re¬ 
ported)  as  readings  once  an  hour  rather  than  as  the  hourly 
mean  that  Adolf  Schmidt  advocated  in  1905  and  that  was 
widely  and  rapidly  adopted  Svalgaard  and  Cliver  [2005] 
showed  that  although  the  variance  of  single  values  is 
larger  than  for  averages,  the  overall  effect  on  IDV  was  small 
(at  most  a  few  percent).  ( lTiis  effect  is  significant  for  the 


HIV  index,  but  for  that  case,  correction  of  the  effect  is 
straightforward  [Sva /guard  and  Cliver ,  2007b).)  The  two 
long-running  series  POT-SFD-NGK  and  PSM-VLJ  CLF 
afford  a  convenient  additional  test  of  this:  POT  changed 
from  values  to  means  with  the  1905  yearbook,  but  CLF 
changed  much  later,  with  the  1972  yearbook,  so  we  can 
directly  compare  the  (raw  uncorrected  in  any  way)  /De¬ 
values  for  the  two  series  (Figure  2).  It  is  evident  that  the 
change  from  hourly  instantaneous  values  to  hourly  means 
did  not  introduce  any  sudden  changes  in  IDV  at  the  times  of 
the  transitions.  The  Japanese  station  at  KAK  changed  from 
values  to  means  in  1955.  The  ratio  between  raw  IDV  for 
KAK  and  SED-NGK  (crosses  on  Figure  2)  also  does  not 
show  any  change  in  1955.  The  American  stations  CLH  and 
I  ION  changed  to  means  with  the  1915  yearbook.  Compar¬ 
ison  (Figure  3)  over  a  24-year  interval  centered  on  1915 
with  the  stations  VLJ  and  DBN,  which  did  not  change 
sampling  procedure,  also  shows  no  detectable  change  in 
IDV  due  to  the  change  in  sampling:  the  ratio  between 
average  CLH- 1  ION  and  VLJ-DBN  is  I  0792  before  1915 
and  1.0792  after  the  change  to  hourly  means  in  1915.  We 
conclude  that  changes  are  too  small  to  justify  attempting  ad 
hoc  correction  based  on  extrapolation  of  modem  data. 

2.1.3.  Using  the  w-Measure  Before  1872 

[n]  Bartels  [1932;  Joos  et  al.,  1952]  compiled  the  it 
measure  from  the  interdiumal  variability  of  the  Horizontal 
Component,  //,  from  hourly  or  daily  values  from  several 
observatories  operating  from  1 872  onwards  as  described  in 
his  paper.  He  wrote,  "Before  1872,  no  satisfactory'  data  for 
the  calculation  of  interdiumal  variabilities  are  available,"  but 
"more  for  illustration  than  for  actual  use."  he  attempted  to 
extend  the  series  backward  to  1835  For  this  he  used  the 
"Einheitliche  Deklinations  Vanationen"  (Unified  Declina¬ 
tion  Variations),  £,  of  Wolf[  1884]  and  the  "summed  ran¬ 
ges,"  .v,  derived  from  the  mean  diurnal  variation  of  //  at 
Colaba  (Bombay)  due  to  Woos  [1910].  lie  derived  regres¬ 
sion  formulae  relating  E  and  s  to  //  for  times  after  1872  and 
used  them  to  synthesize  values  of  u  for  the  earlier  years; 
giving  v  double  the  weight  of  E.  Bartels  justified  this  by 
showing  that  for  the  annual  means  1872-1901,  the  values  of 
//  derived  from  II  and  the  values  of  .v  have  a  high  linear 
correlation  coefficient.  (From  E  and  we  calculate  a  value 
of  0.72  for  the  value  for  it  for  1857.5  using  the  formulae 
given  by  Bartels.)  We  have  extended  his  analysis  by  cal¬ 
culating  the  correlation  between  IDV  and  the  Summed 
Ranges  for  1872-1905  finding  a  correlation  coefficient  of 
0  86.  Figure  4  shows  the  agreement  between  observed  IDV 
(red)  and  that  calculated  from  s  (blue). 

[12]  Furthermore,  as  shown  by  Svalgaard  and  Cliver 
[2005]  there  is  a  good  linear  correlation  between  IDV  (or 
1IMF  B  derived  from  it)  and  the  square  root  of  the  sunspot 
number,  R  The  main  sources  of  the  equatorial  components 
of  the  Sun’s  large-scale  magnetic  field  are  large  active 
regions.  If  these  active  regions  emerge  at  random  longitudes, 
their  net  equatorial  dipole  moment  will  scale  as  the  square 
root  of  their  number.  Thus  their  contribution  to  the  average 
I  IMF  strength  will  tend  to  increase  as  Rl  “  (for  a  detailed 
discussion,  see  Wang  and  Sit ee lev  [2003]  and  l Vang  et  al. 
[2005]). 

[13]  To  the  extent  that  the  //-measure  before  1872  can  be 
taken  as  a  geomagnetic -based  measure  of  the  sunspot 
number,  it  is  therefore  to  be  expected  that  the  //-measure 
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No  Significant  Effect  Due  to  Changing  from  Hourly  Values  to  Hourly  Means 


Figure  2.  IDV  calculated  without  any  normalization  or  adjustments  for  the  long-running  German  series 
(Potsdam  POT-Seddin  SED-Niemegk  NGK;  red  curves)  and  the  long-running  French  series  (Parc  Saint- 
Maur  PSM-Val  Joyeux  VLJ-Chambon-la-Foret  CLF;  green  curves).  Vertical  lines  show  when  the 
replacement  stations  went  into  operation  and  the  ovals  show  when  the  yearbook  values  changed  from 
being  instantaneous  hourly  spot  values  to  hourly  means.  The  blue  (spot  values)  and  red  (hourly  means) 
crosses  show  the  ratio  between  raw  IDV s  for  KAK-Kakioka  and  SED-NGK.  KAK  changed  from  record¬ 
ing  spot  values  to  recording  hourly  means  in  1955.  There  are  no  clear  indications  of  changes  in  IDV  due  to 
the  change  in  recording  reporting  practice. 


also  will  serve  as  a  proxy  for  IDV.  This  estimate  will  be 
independent  of  any  assumptions  about  the  constancy  of  the 
calibration  of  the  sunspot  number  (compare  the  difference 
between  the  Zurich  Sunspot  Number  and  the  Group  Sunspot 
Number  [Hoyt  et  al. ,  1994]) 

[u]  Figure  5  shows  that  /DU can  also  be  directly  inferred 
from  the  daily  range,  rK,  of  the  east  component  (equivalent 
to  the  Declination  for  this  purpose)  of  the  geomagnetic  field 
and  that  therefore,  again,  that  the  //-index  before  1872 
(strongly  influenced  by  the  range  of  the  daily  variation)  can 
be  used  for  estimation  of  IDV ,  albeit  with  slightly  less 
accuracy  than  after  1872.  This  conclusion  may  seem  at 
variance  (and  did  surprise  us)  with  our  initial  decision  to  use 
only  nighttime  data  in  the  derivation  of  IDV ,  but  emerges 


naturally  (and  inescapably)  after  our  analysis  had  shown  that 
IDV  derived  without  any  dependence  on  daytime  data  is 
comparable  to  IDV  derived  from  daily  ranges  because  of  the 
strong  dependence  of  both  on  the  sunspot  number.  This  is 
clearly  demonstrated  in  Figure  6  that  shows  raw  IDV  cal¬ 
culated  for  PSM-VLJ-CLF  and  POT-SFD-NGK  deter¬ 
mined  from  nighttime  differences  (blue)  and  daytime 
differences  (red).  This  realization  opens  the  door  for  use  of 
19th  geomagnetic  stations  that  only  observed  during  the  day 
as  long  as  the  observations  were  made  at  fixed  hours. 

[15]  For  the  reasons  given  above,  we  find  that  IDV  can  be 
estimated  with  confidence  from  Bartels’  //-measure  also 
before  1872,  justify  ing  our  reconstruction  of  I  IMF  B  since 


No  Difference  Between  Spot  and  Mean  Values  Before  and  After  1915 


Figure  3.  Raw  IDV  for  the  average  of  stations  VLJ  and  DBN  (both  reporting  instantaneous  “spot"  values 
ev  ery  hour  for  the  12-year  intervals  before  1915  and  after  1915)  and  for  the  average  of  CLH  and  I  ION 
(reporting  spot  values  before  1915  (pink)  and  hourly  mean  value  thereafter  (blue)) 
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Comparison  of  IDV  and  Colaba  Summed  Ranges 


Figure  4.  Comparison  of  observed  IDV{  red  open  squares)  and  synthetic  IDV  calculated  from  s  using  the 
regression  equation  shown  (solid  blue  circles)  derived  from  the  correlation  between  IDV  and  the  Summed 
Ranges,  v,  of//  from  Colaba  [Moos.  1910;  p.  294,  Table  261]  for  1872-  1905. 


1835,  with  expected  refinement  of  the  reconstruction  upon 
further  digitization  of  19th  century  geomagnetic  yearbooks. 
2.1.4.  /DC- Index  1835-2009 
[i6]  From  the  -1,375,000  daily  differences  (3775  station- 
years)  denved  from  the  stations  in  Table  1  we  construct  the 
/DC-index  shown  in  Figure  7,  with  individual  station  curves 
in  gray.  The  composite  (red)  curve  is  the  mean  of  the 
median  and  average  values  for  each  year,  while  before  1872 
the  dashed  curve  shows  IDV  estimated  from  //.  Also  shown 
(blue  curve)  is  the  number  of  stations  contributing  to  the 
mean.  The  large  number  of  stations  from  1957  on  does  not 


add  further  significance  to  the  composite,  but  only  serves  to 
establish  the  range  of  scatter  of  the  values. 

[17]  It  is  evident  that  IDV  from  only  a  single  station 
(provided  that  not  too  much  data  is  missing  either  because 
the  recording  went  off  scale  or  as  a  result  of  other  problems) 
does  not  differ  much  from  the  mean  of  many  stations;  the 
standard  deviation  of  /DT-values  for  all  stations  for  a  given 
year  is  less  than  I  nT  or  about  9%.  This  means  that  only  a 
few  (good)  stations  are  needed  for  a  robust  determ i nation  of 
IDV.  This  conclusion,  of  course,  only  emerges  after  the 
spread  of  /DT- values  has  first  been  shown  to  be  small.  The 


IDV  as  a  Function  of  Daily  Range  of  East  Component 


Figure  5.  IDV  plotted  against  the  amplitude  of  the  daily  range,  rY ’,  of  the  east  component  (Data  Set  S2  in 
the  auxiliary  material)  of  the  geomagnetic  field  denved  from  PSM-VLJ-CLF  and  POT-SED-NGK,  cov¬ 
ering  the  interval  1883  2008  (dark  blue  diamonds)  for  which  we  have  data  for  these  stations.  Since  1957. 
the  number  of  stations  contributing  to  IDl  is  high  ( —50)  for  every  year,  so  the  data  is  good.  The  open  red 
squares  show  the  same  relationship  for  1957  2008. 
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IDV  Derived  Using  Nighttime  (blue)  and  Daytime  (red)  Differences 


Figure  6.  Raw  IDV  derived  using  nighttime  hourly  data  (blue)  and  daytime  hourly  data  (red)  for  the 
French  stations  PSM-VLJ-CLF  and  the  German  stations  POT-SED-NGK.  The  daytime  values  are 
~~30%  larger  than  the  nighttime  v  alue  because  of  day-to-day  (non-solar)  variations  of  the  regular  solar 
variations,  SR%  and  have  been  scaled  to  the  same  average  as  the  nighttime  values  for  easier  comparison. 


standard  error  of  the  mean  of  more  than  fifty  stations  is 
0.1  nT. 

[is]  Figure  8  shows  that  the  differences  between  1DV05 
and  1DV09  are  slight,  and  due  to  the  additional  data  since 
1880.  During  the  period  of  overlap  (1872-2003,  2004  was 
only  partial),  the  two  time  series  agree  within  an  RMS  of 
0.33  nT  or  3%.  The  coefficient  of  determination  for  the 
correlation  between  IDV09  and  1DV05  is  R2  =  0.984.  IDV  is 
a  robust  index. 


2.1.5.  Physical  Interpretation  of  IDV:  Measure  of  the 
Energy  in  the  Ring  Current 

[w]  Svalgaard  and  Cliver  [2005]  reported  that  IDV  is 
closely  correlated  with  the  negative  part  of  the  Ds(- index 
based  on  data  back  to  1932  [Karmen  and  Mursula,  2005]. 
Svalgaard  and  Cliver  [2006]  extended  that  relationship  back 
to  1905  using  the  100-year  Dv/-series  derived  by  Love 
[2007;  J.  J.  Love,  personal  communication,  2006],  and 
confirm  it  here  using  IDV09.  Yearly  averages  of  Dv,  (scaled 
to  Kyoto  Dst\  we  use  Dv,  here  in  a  generic  sense  without 


Figure  7.  Yearly  IDV-  indices  derived  for  individual  stations  (as  given  in  Table  1 )  shown  as  gray  curves. 
I  Tie  red  curve  is  a  composite  index  calculated  as  the  mean  of  the  median  and  average  values  of  the 
individual  station  values.  This  procedure  may  be  justified  by  the  very  small  difference  between  medians 
and  averages  (0. 16  nT  on  av  erage,  see  Figure  8).  The  number,  N%  of  contributing  stations  is  show  n  by  the 
thin  blue  curve  and  the  corresponding  number  for  IDV05  as  a  dashed  light  blue  curve  The  //-measure  is 
considered  a  single  station.  A  few  station  values  differing  more  than  five  standard  deviations  from  the 
average  for  a  given  year  were  omitted  in  calculating  the  average  for  that  year,  fhe  dashed  tune  shows 
IDV  derived  from  the  //  measure. 
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Details  of  Evolution  of  IDV 


1830  1840  1850  1860  1870  1880  1890  1900  1910  1920  1930  1940  1950  1960  1970  1980  1990  2000  2010 

Figure  8.  Average  yearly  values  of  IDV09  (dark  blue  curve)  compared  with  median  yearly  values  (light 
blue  curve)  and  compared  with  published  1DV05  (red  curve). 


distinguishing  between  different  derivations  of  the  under¬ 
lying  physical  measure  sought  captured  by  Dst)  when  the 
hourly  value  was  negative  were  computed  and  found  to  be 
strongly  correlated  with  IDV  { R2  -0.91):  IDV=  -0.45  (D„< 
0).  Figure  9  compares  1DV09  and  IDV  computed  from  Dsl. 
The  good  match  suggests  that  IDV  is  a  measure  of  the  same 
physical  reality  as  negative  Dv/,  namely  the  energy  in  the  Ring 
Current,  which  then  in  turn  seems  to  be  controlled  by  HMF  B. 
{Dsl  <0)  4.81  B  -9.41  (R2  =  0.84),  and  we  can  then  also  use 
Dxt  to  determine  the  l  IMF  strength:  B  —  2.70  -  0. 1 736  (Dst  < 
0).  Schmidt  [1926]  actually  suggested  that  in  the  definition  of 
the  //  measure  it  would  be  slightly  better  to  only  use  the 
negative  differences  between  consecutive  days. 

2.2.  Using  IDV09  to  Calculate  HMF  Strength, 
1835-2009 

[20]  Since  the  2005  definition  paper,  lower  values  of 
I  IMF  strength,  /?,  have  improved  the  dy  namic  range  (and 


thus  the  statistical  significance)  of  the  correlation  between 
IDV'  and  B  An  approximate  linear  correlation  was  found, 
but  there  is  no  a  prion  reason  the  relationship  would  be 
strictly  linear  In  addition,  it  has  been  argued  [Lockwood 
et  al .,  2006]  that  B  should  be  taken  as  the  independent 
variable  instead  of  IDV  Svalgaard  and  Cliver  [2006]  showed 
that  it  does  not  make  much  difference  which  way  the  corre¬ 
lation  is  evaluated.  In  the  end,  the  RMS  difference  (0.4  nT  or 
less  than  -10%)  between  HMF  B  observed  in  situ  near  the 
Earth  and  inferred  from  IDV'  is  what  matters.  (Using  hourly 
av  erages  from  the  OMNI  dataset  for  histoncal  data  and  from 
ACE  for  near  real  time  recent  data.)  The  average  coefficients 
for  the  linear  correlation  performed  four  ways  (average, 
median,  and  for  each  direct  and  inverse)  are 

fl(nT)  =  (2.06  ±0.21)  +  (0  441  ±  0.021  )IDV 

(R2  =  0.869)  (2) 


Negative  D*t  and  IDV  measure  the  same  effect 


Figure  9.  Yearly  average  values  of  IDV  mA  of  D„  when  it  was  less  than  zero  (based  on  Dsl  from  Kyoto 
W  DC  and  on  Dst  from  J.  J.  Love  (personal  communication,  2006)  scaled  to  Kyoto  levels).  The  “spike”  in 
1 909  is  due  to  the  extremely  strong  storm  on  25  September  1909  causing  loss  of  data  at  all  but  one  station 
(  \PI),  giv  ing  that  one  data  point  undue  influence.  To  guard  against  the  influence  of  such  sporadic  extreme 
values,  ihe  daily  values  of  IDV  were  capped  at  75  nT. 
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Table  2.  IDV- Index  Measured  or  Inferred  for  Each  Year  Since  Table  2.  (continued) 


Year 

1835.5 

1836.5 

1837.5 

1838.5 

1839.5 

1 840.5 

1841.5 

1842.5 

1843.5 

1844  5 

1845  5 

1846  5 

1847.5 
1848  5 

1849.5 
1 850  5 

1851.5 

1852.5 

1853.5 

1854.5 

1855.5 

1856.5 

1857.5 

1858.5 

1859.5 

1860.5 

1861.5 

1862.5 

1863.5 

1864.5 

1865.5 

1866.5 

1867.5 

1868.5 

1869.5 
1870  5 

1871.5 

1872.5 

1 873.5 

1874.5 

1875.5 

1876.5 

1877.5 

1878.5 

1879.5 
1880  5 

1881.5 

1882.5 
1883  5 

1884.5 

1885.5 

1886.5 

1887.5 

1888.5 

1889.5 
I  890.5 
189!  5 
1892  5 

1893.5 
I  894.5 

1895.5 

1896.5 

1897.5 

1898.5 
1899  5 

1900.5 

1 90 1. 5 

1902.5 
1903  5 


IDV  09 


11.60 
16.30 
I  6.00 
16.80 
14.00 
12.20 
10.10 
9.00 

8.90 
8.50 
9  50 
10.10 

1 1.40 

13.10 

12.00 

9.90 
9  00 
7.60 

7.80 

7.60 

5.80 

6.60 

7.20 
10.60 

14.30 
13.50 

12.20 
10.00 

9.40 

8.40 
7  90 

7.30 

6.80 

8.10 
11.10 
16  70 
1 6.00 

14.60 
9  90 
9  50 

7.20 
6  00 

6.60 
5.80 
6  00 
7.89 
8.69 
12.47 

10.20 
9.36 
9.57 
9.14 

7.75 
727 
6.99 
6.22 
8  60 
14  02 
10.79 
13.12 
9  95 
10.48 
8.71 
8.98 
7  06 

5.75 
4.90 
5.04 
7  03 


Observed  Vear 


IDV  HMF  B 

HMF  B 

7.23 

I  VU4  S 

1905  5 

9.30 

1906.5 

9  17 

1907.5 

9.51 

1908  5 

8  30 

1909.5 

7.50 

1910.5 

6.55 

191  1.5 

6.05 

1912.5 

6  00 

1913.5 

5.81 

1914.5 

6.28 

1915.5 

6.55 

1916.5 

7.14 

1917.5 

7.90 

1918.5 

7.41 

1919.5 

6.46 

1920.5 

6  05 

1921.5 

5  39 

1922  5 

5  48 

1923.5 

5  39 

1924.5 

4.51 

1925.5 

4.90 

1926  5 

5  19 

1927.5 

6.78 

1928.5 

8.43 

1929  5 

808 

19305 

7.50 

1931.5 

6  51 

1932.5 

6.23 

1933.5 

5.76 

1934.5 

5.53 

1935.5 

5.24 

1936.5 

5.00 

1937.5 

562 

1938.5 

7.01 

1939.5 

9  47 

1940.5 

9  17 

1941.5 

8.56 

1942.5 

6  46 

1943.5 

6.28 

1944  5 

5  19 

1945  5 

4  61 

1946  5 

490 

1947.5 

4.51 

1948  5 

4.61 

1949.5 

5.53 

1950.5 

5.90 

1951  5 

7.62 

1952  5 

6.60 

1953.5 

6.22 

19545 

6.22 

1955.5 

6.11 

1956.5 

5.46 

1957  5 

5  23 

1958.5 

5  09 

1959.5 

4.72 

1960  5 

5.86 

1961  5 

8  31 

1962  5 

6.87 

1963.5 

791 

1964.5 

6  48 

1965.5 

6.73 

1966.5 

591 

1967.5 

6.04 

1968.5 

5.13 

19695 

4.49 

1970  5 

4.06 

1971  5 

4  13 

1972.5 

5  11 

1971  5 
1974.5 

IDV  09 


7.54 
8.62 
7.49 
8  83 
9.45 
9.95 
8.10 
7.08 
5.69 
5.15 
6.22 
8.09 
9.19 
10  95 

10.97 
11.57 
10.28 

8.97 

7.74 
6.17 
6.89 
8.05 

10.69 
9  29 

9.69 
9.64 
10.22 
7.38 
7  22 
6  96 
6.83 

7.75 
8.81 
12.1  I 
14.02 

12  79 
12.48 
12.10 
9.57 
8.97 
8.28 

8.75 
14  33 
13.85 
10.87 

13  55 
12  56 
12.46 
10.97 
8.90 
746 
8.69 

13.38 
16.65 
15.42 

14.39 
15.87 
1  1.49 
8.62 
8.06 
7.19 
6.93 
7.88 
10.30 
9.47 

9.39 
10.12 
8.84 
9  49 
9.28 

9  18 


HW  HMF  B 


5.36 

5.87 

5.33 

5  97 

6  26 
6.48 

5.63 
5.14 
4.46 

4.18 
4.72 
5.62 
6  13 

6.94 

6.95 
7.22 

6.64 
6.03 

5.45 
4.70 
5.04 

5.60 
6.82 

6.18 
6.37 
6.34 

6.61 
5.28 
5.21 
5.08 
5.02 

5.46 
5.96 

7.46 

8.31 
7,77 
7.63 
7.46 

6.31 
603 
5.71 
5.93 
8  44 
8  24 

6  91 
8  10 

7  66 
762 
6.95 
6.00 
5  32 
5.90 
802 
9.45 
8.92 
8.47 
9.11 
7  18 
5.87 
5.60 

5  19 
5.07 
5.52 

6  65 
6  26 
6.23 
6  56 

5  97 

6  27 
6  18 
6. 1  3 


Observed 
HMF  B 


5  45 

5  12 
5.06 
6.00 

6  36 
6  19 
6  05 
6.35 
6.00 
6.38 
6  35 
6.63 
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Table  2.  (continued) 


Year 

IDV  09 

nn  HMF  B 

Observed 

HMF  B 

1975.5 

8.15 

5.65 

5,82 

1976.5 

8.70 

5  91 

5.45 

1977.5 

8.96 

6.03 

5.85 

1  978.5 

12  32 

7.56 

7  08 

1979.5 

11.78 

732 

7.59 

1980.5 

10.51 

6.74 

6  98 

1981.5 

13.78 

8  20 

7.84 

1982.5 

15.25 

8.84 

881 

1981  5 

11.60 

7.23 

7.61 

1984.5 

10.50 

6.74 

7.52 

1985.5 

9.06 

6.07 

5.89 

1986.5 

8.80 

5.95 

5.74 

1987.5 

8.20 

5.67 

6.09 

1988.5 

10  21 

6.61 

7.30 

1989.5 

1674 

9.48 

8.15 

1990.5 

12  84 

7  79 

7  29 

1991.5 

15.77 

9.07 

9.34 

1992.5 

12  87 

7.80 

8.25 

199}.  5 

10.08 

6.54 

6.59 

1994.5 

9  06 

6.07 

6.15 

1995.5 

9.08 

6  08 

5.72 

1996  5 

6.76 

4.98 

5.11 

1997.5 

8.06 

5.60 

5.51 

1998  5 

10.34 

6.66 

6.89 

1999  5 

9.82 

6.42 

6.91 

2000.5 

13  36 

8.02 

7.18 

2001.5 

13  44 

8.05 

6.94 

2002.5 

10.90 

6.92 

7.64 

2003.5 

12  51 

764 

7  60 

2004.5 

9.42 

6  24 

6.53 

2005.5 

9  44 

6.25 

6.25 

2006  5 

7.22 

5.21 

5.03 

2007  5 

5  96 

4  59 

4  48 

2008  5 

5  29 

4  25 

4.23 

2009.5 

5.04 

4.13 

4.05 

2010  2 

5  50 

4.45 

4.95 

‘The  HMF  strength  flat  the  Earth  is  derived  from  IDV as  per  section  2.2. 
The  field  observed  in  situ  (OMNI  and  ACE  data  sets)  is  given  for  com¬ 
parison  A  few  years  had  very  incomplete  data  coverage  and  missing  data 
were  derived  by  linear  mlerpolation  across  data  gaps  to  avoid  uneven 
coverage  skewing  Ihe  average  Those  \alucs  are  in  italics. 


The  equivalent  power  law  dependence  comes  to 

5(nT)  (1.33  ±0.07)/DFH,WWtM02,(R:  0.905)  (3) 

The  adopted  values  for  B  inferred  from  1DV09  given  in 
Table  2  are  the  mean  values  calculated  using  these  two 
relationships.  Data  Set  S3  in  the  auxiliary  material  sum¬ 
marizes  the  coefficients  for  all  correlations.  The  ‘error  bars' 
quoted  are  not  a  measure  of  the  statistical  significance  of  the 
correlations  in  a  strict  sense,  but  are  solely  indicative  of  the 
range  or  variability  of  the  various  coefficients. 

[21 J  Figure  10  shows  the  values  for  HMF  B  inferred  from 
IDV  from  1835  to  the  present  (blue  curve)  and  B  measured 
by  spacecraft  (red  curve).  A  4th  order  polynomial  tit  sug¬ 
gests  a  -100  year  Gleissberg  cycle.  Cycle  23  looks 
remarkably  like  cycle  13,  including  the  very  deep  solar 
minimum  following  both  cycles,  likely  presaging  a  weak 
cycle  24  as  predicted  from  the  solar  polar  fields  [Svalgaard 
et  al .,  2005;  Schatten ,  2005].  It  is  clear  that  we  are  returning 
to  conditions  prevailing  a  century  ago.  It  seems  likely  that 
other  solar  parameters  such  as  Total  Solar  Irradiance 
[Frohlich,  2009;  Steinhilber  et  al.,  2009]  and  cosmic  ray 
modulation  [Steinhilber  et  al .,  2010]  are  reverting  to  similar 
conditions  with  possible  implications  for  the  climate-change 
debate. 

2.3.  Comparison  of  IDV09-Bascd  B  With  Other  Recent 
Reconstructions 

2.3.1.  Consilience  of  Reconstructions  Based  on 
Geomagnetic  Data 

[22]  Reconstructions  of  1IMF  B  have  been  discordant  in 
the  past  (e.g.,  Svalgaard  [1977],  as  referenced  by  Schatten 
et  al  [1978],  Andreasen  [1997],  Lockwood  et  al  [1999, 
2006],  and  Svalgaard  and  Cliver  [2005,  2006,  2007b)).  The 
realization  [Svalgaard  et  al.,  2003]  that  geomagnetic  indices 
can  be  constructed  that  have  different  dependencies  on  B 
and  solar  wind  speed  (F)  has  enabled  robust  determinations 
of  both  V  [Svalgaard  and  Cliver ,  2007b;  Rouillard  et  al., 
2007;  Lockwood  et  al . ,  2009]  and  B  [Svalgaard  and 
Cliver ,  2005,  2006;  Lockwood  et  al ,  2009]  that  have  con¬ 
verged  to  a  common,  well -constrained  data  set.  Progress  has 


B  nT  Heliospheric  Magnetic  Field  Strength  B  (at  Earth)  Inferred  from  IDV  and  Observed 


Year 


Figure  10.  Yearly  average  values  of  the  I  IMF  B  inferred  from  the  IDV  index  (dark  blue  curve)  and  from 
1  he  early  a  measure  (light  blue  curve)  compared  with  in  situ  measurements  (red  curve).  There  is  a  hint  of 
ihe  -100  year  Gleissberg  cycle. 
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B  nT  Heliospheric  Magnetic  Field  Strength  B  (at  Earth)  Inferred  from  IDV  and  Observed 


Figure  1 1.  Comparison  of  I  IMF  B  determined  from  IDV  (light  blue  curve  using  equation  (3)  and  dark 
blue  curve  using  equation  (2)),  by  Lockwood  et  ai.  [2009]  (green  curve),  and  observed  by  spacecraft  (red 
curve). 


been  swift  and  Figure  I  I  shows  the  convergence  of  HMF  B 
determined  by  Lockwood  et  al.  [2009]  to  the  values  deter¬ 
mined  from  IDV  [Svalgaard  and  C liver ,  2005;  this  paper]. 
The  Lockwood  ct  al.  [2009,  and  references  therein]  recon¬ 
struction  is  still  discordant  from  ours  for  a  few  years  during 


solar  cycle  14,  but  apart  from  that,  the  agreement  is  quite 
remarkable  and  the  issues  seem  resolved. 

[21]  Figure  12  details  the  evolution  of  the  various 
determinations  of  B  since  the  influential,  but  now  superseded, 
Lockwood  et  al.  [1999]  paper.  It  is  clear  that  we  now  possess 
the  methodology  to  infer  B  with  good  accuracy  as  far  back 


Convergence  of  HMF  B  of  Lockwood  &  Rouillard  et  al.  to  Svalgaard  &  Cliver 


Figure  12.  Comparison  between  I  IMF  B  derived  by  Svalgaard  and  Cliver  [2005]  (light  blue  curve  and 
open  circles);  this  paper  (dark  blue  curve  and  open  circles);  and  I  IMF  B  derived  by  Lockwood  et  al. 
[1999]  (orange  curve  and  pluses),  Rouillard  et  al  [2007]  (their  point  for  1901  was  in  error)  and  A. 
Rouillard  (personal  communication,  2007)  (pink  curve  and  pluses),  and  Lockwood  et  al.  [2009)  (red  curve 
and  pluses),  matched  to  in  situ  observations  of  B  (black  dots). 
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Heliospheric  Magnetic  F  eld  Comparisons 


Figure  13.  Yearly  averages  of  near-Farth  11MF  B  inferred  by  this  paper  (blue  curve  /?s&t ).  by 
Lockwood  et  ai  [2009]  (green  curve  tfrR&r),  observed  by  spacecraft  (red  curve  Z?0bs)  compared  to  B 
inferred  by  McCracken  [2007]  (purple  curve  Z?McC)  The  large  arrow  marks  the  beginning  of  the  neutron 
monitor-based  part  of  the  record  One  might  speculate  that  the  extremely  low  values  during  1883-1893 
are  caused  by  the  explosion  of  Krakatoa  ejecting  sulfur-rich  aerosols  into  the  stratosphere  influencing  the 
deposition  of  H  Be. 


as  continuous  geomagnetic  records  of  H  reach.  A  concerted 
effort  of  digitization  of  19th  century  yearbook  records  would 
promise  to  further  improve  our  know  ledge  of  the  magnetic 
field  in  the  heliosphere 

[24]  Svalgaard  and  Cliver  [2007a]  argued  for  a  floor  in 
yearly  averages  of  solar  wind  B  which  was  approached  at 
every  1 1  yr  minimum  and  represented  the  ground-state  of 
the  Sun  during  extended  minima  such  as  the  Maunder 
Minimum.  With  the  larger  dynamic  range  afforded  by  the 
current  minimum,  we  can  now  refine  the  value  of  the  floor 
to  be  closer  to  the  ~4  nT  observed  during  2008  and  2009 
[see  also  Owens  et  al. ,  2008],  returning  to  the  values 
inferred  for  11 -yr  minima  during  the  previous  Gleissberg 
minimum  at  the  turn  of  the  20th  century. 

2.3.2.  Comparison  With  l0Be- Based  Reconstructions 

[25]  McCracken  [2007]  spliced  together  I0Be  data,  ioni¬ 
zation-chamber  cosmic  ray  data  (calibrated  with  balloon 
flight  data),  and  neutron  monitor  cosmic  ray  data  to  produce 
an  ‘equivalent’  neutron  monitor  count  series  covering  the 
entire  interval  1428-2005,  and  inverted  the  senes  for  B  in 


order  to  express  the  data  in  terms  of  the  1  IMF  B.  In  Figure  1 3 
we  compare  his  series  for  HMF  B  with  the  ‘consensus’  B 
from  geomagnetic  data. 

[26]  In  McCracken’s  time  series  for  B ,  a  large  step- like 
change  (1.7  nT;  from  3.5  nT  to  5.2  nT;  the  largest  jump  in 
the  entire  -600-year  record)  occurs  between  the  1944  and 
1954  sunspot  minima  flanking  cycle  18.  No  such 
corresponding  change  is  observed  in  the  concordant 
reconstructions  of  Svalgaard  and  Cliver  [2005;  this  paper], 
Ronillard  et  ai  [2007],  and  Lockwood  et  al.  [2009],  nor  in  B 
calculated  from  the  quantity  fll'deduced  by  Le  Sager  and 
Svalgaard  [2004]  using  either  Fof  Svalgaard  and  Cliver 
[2006]  or  of  Ronillard  et  al  [2007],  or  in  B  deduced  from  Dst. 

[27]  Mnscheler  et  al.  [2007]  discuss  the  uncertainties  w  ith 
the  balloon-borne  data  that  form  the  basis  for  McCracken’s 
calibration  of  the  composite  equivalent  neutron  monitor  data 
before  1951.  The  strong  geomagnetic  evidence  argues  that 
this  particular  calibration  of  the  pre-neutron  monitor  cosmic 
ray  reconstruction  is  not  on  a  firm  footing.  We  suggest  that 
part  of  the  reason  for  the  disagreement  might  lie  with  the 


Svalgaard 


McCracken(2607) 

&  Cllv< 
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--  a=1  a=2  .  — <^Avarying  vsw±u 

1750  1800  1850  1900  1950  201 


Year 


00° 


Figure  14.  Comparison  of  our  reconstruction  of  1  IMF  B  (red  curve,  25  year  running  means)  with  other 
reconstructions  as  reported  by  Steinhilher  et  al.  [2010.  Figure  7)  (adapted  and  reproduced  here),  e.g.,  w  ith 
their  25-year  running  mean  of  their  PC  A  based  reconstruction  (purple  dashed  curve)  of  B.  The  oval 
outlines  an  area  of  disagreement  for  which  sufficient  geomagnetic  data  exists  that  may  be  used  to  resolve 
rhe  discrepancy 
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IDV  and  HMF  B  1880-1920  (All  stations) 


Figure  15.  IDV  (blue  curves)  and  inferred  HMF  B  (red  curve;  dashed  line:  25-year  running  mean)  1 880- 
1020  for  all  stations  (as  noted  by  their  1AGA  designations  -  10//  shown  as  a  dashed  gray  curve)  where 
good  geomagnetic  data  are  available  so  far. 


calibration  and  splicing  together  of  the  disparate  cosmic  ray 
data  sets. 

[:x]  After  our  paper  was  submitted,  we  were  pleased  to 
read  a  paper  by  Steinhilher  et  al.  [2010]  in  which  a  new 
"Be-based  reconstruction  has  moved  closer  to  our  recon¬ 
struction,  to  that  of  Rouillard  et  al.  [2007],  and  to  that  of 
Caballero -Lopez  et  al.  [2004]  w  ith  diffusion  coefficient 
depending  inversely  of  B~  (a  -  2).  The  reconstruction  of 
Steinhilher  et  al.  [2010]  still  differs  somewhat  with  the 
geomagnetic  based  reconstructions,  especially  for  the 
-1880-1900  interval  (Figure  14)  and,  just  like  the  previous 
discrepancy,  this  will  need  to  be  resolved.  We  suggest  that  if 
the  sharp  dip  around  -1895  is  not  borne  out  by  further 
investigation,  the  magnitude  of  earlier  excursions  to  very 
low  values  may  also  be  in  doubt.  Figure  15  shows  IDV  for 
all  stations  for  the  interval  1880-1920  and  does  not  support 
the  marked  decrease  around  -1895.  It  is  unlikely  that  data 
from  further  stations  will  change  that  conclusion.  In  line 
with  this  conclusion,  W.  R.  Webber  et  al.  (A  comparison  of 
new  calculations  of  the  yearly  ,0Be  production  in  the  Earth’s 
polar  atmosphere  by  cosmic  rays  with  yearly  uBe  mea¬ 
surements  in  multiple  Greenland  ice  cores  between  1939 
and  1994— A  troubling  lack  of  concordance,  manuscript  in 
preparation,  2010)  suggest  that  “more  than  50%  of  the  nBe 
flux  increase  around,  e.g.,  1700  A.D.,  1810  A.D.  and  1895 
A.D.  is  due  to  nonproduction  related  increases." 

3.  Summary  and  Discussion 

[:*>]  We  have  extended  our  1872  2004  11MF  time  series 
[Svalgaard  and  Cliver ,  2005]  to  the  years  1835-2009 
(Figure  10).  The  1835-1871  interval  is  based  on  Bartels’ 
//  measure,  which  he  extended  from  1871  back  to  1835 
using  Wolfs  [1884]  Declination  index  based  on  several 
European  stations  and  Moos'  [1910]  Summed  Ranges  from 
Colaba.  The  1872  2009  interval  is  based  on  the  IDV- index, 
with  significantly  more  data  for  the  early  years  (1872- 
1910).  The  forward  extension  of  the  11MF  series  through 
2009  is  important  because  the  years  2007  2009  witnessed 


the  lowest  annual  averages  of  IDV  during  the  space  age.  For 
the  time  of  overlap  between  the  re-evaluated  IDV- index 
(1DV09)  and  1DV05,  the  difference  is  very  small,  testifying 
to  the  robustness  of  the  index. 

[30]  A  comparison  of  lDV09-based  HMF  strength  with 
those  obtained  by  other  investigators  using  various  combi¬ 
nations  and  permutations  of  geomagnetic  indices  revealed  a 
pleasing  agreement  (Figure  1 1)  in  what  had  been  previously 
a  contentious  field  of  research  (Figures  12  and  13).  The 
technique  proposed  by  Svalgaard  et  al.  [2003]  and  adopted 
by  Rouillard  et  al.  [2007]  to  use  indices  with  different 
dependencies  on  B  and  V  to  separate  these  variables  has 
proven  out  and  allowed  the  vast  storehouse  of  hourly  and 
daily  data  to  be  brought  to  bear.  In  particular,  the  B  values 
deduced  and  cross-checked  [Le  Sager  and  Svalgaard,  2004] 
by  this  method  have  substantiated  the  approach  made 
possible  by  the  IDV-  index  as  well  as,  as  Svalgaard  and  Cliver 
[2005]  suggested  and  has  confirmed  here,  the  negative 
component  of  the  Dvf-  index  (Figure  9).  We  conclude  that  the 
long-term  variation  of  heliospheric  B  is  firmly  constrained 
during  the  time  for  which  it  is  based  on  hourly  values  of  II, 
and  that  current  values  at  the  solar  minimum  between  cycles 
23  and  24  are  back  to  where  they  were  108  years  ago  at  the 
solar  minimum  between  cycles  13  and  14. 

[31]  Although  the  recent  reconstruction  of  B  based  on 

Be  data  [Steinhilher  et  al.,  2010]  generally  agrees  well 

with  the  geomagnetic-based  reconstruction  there  is  dis¬ 
agreement  for  the  decade  just  prior  to  1900  (Figure  14). 
Further  examination  of  these  years  is  critical  because  they 
present  the  only  example  during  the  175  year  interval  of 
geomagnetic-based  B  where  the  floor  [Svalgaard  and 
Cliver ,  2007a]  in  the  solar  wind  is  challenged  by 
Steinhilher  et  al,  [2010],  a  challenge  not  supported  by  the 
geomagnetic  evidence  (Figure  15). 

[32]  Veknowtedgments.  Geomagnetic  data  has  been  downloaded 
from  the  World  Data  Centers  for  Geomagnetism  in  Kyoto.  Japan,  and  C  open- 
hagen.  Denmark  (now  defunct),  and  from  IN  1 1  RM  AGM  T  tl  hop:  www 
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intcniiiignct.org  Data  c.html.  1  he  research  results  presented  in  this  paper 
rely  on  the  data  collected  at  magnetic  observatories  worldwide,  and  we 
thank  the  national  institutions  that  support  them.  We  also  recognize  the  role 
of  the  IN  I  I'.RM  \GNTT  program  in  promoting  high  standards  of  magnetic 
observatory  practice  We  thank  the  many  people  worldwide  who  have 
helped  us  with  collection  of  data  and  metadata  in  addition  to  what  is  avail¬ 
able  from  public  sources.  Wc  thank  Vladimir  Pnpitashvili  for  ihe  program 
to  calculate  corrected  geomagnetic  coordinates  using  the  GUFMI  coeffi¬ 
cients  (courtesy  of  Catherine  Constable).  I  he  OMNI  data  set  was  down¬ 
loaded  from  http  omniweb.gsfc.nasa.gov/.  Real  time  AC  I  inierpla notary 
data  is  downloaded  from  http  wvvw  svvpc  noaa.gov  ftpinena lists  ace2. 
html. 

[.13]  Philippa  Browning  thanks  the  rev  iewers  for  their  assistance  in 
evaluating  this  manuscript. 
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